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Abstract 

We consider phenomenologically allowed structures of the neutrino mass matrix 
in the case of three light neutrino species. Constraints from the solar, atmospheric 
and reactor neutrino experiments as well as those from the neutrinoless double beta 
decay are taken into account. Both hierarchical and quasi-degenerate neutrino mass 
cases are studied. Assuming maximal — i^r mixing we derive simple approximate 
expressions giving the values of the neutrino masses and remaining lepton mixing 
angles in terms of the entries of the neutrino mass matrix. Special attention is paid 
to the small entries which are usually not specified in discussions of the neutrino mass 
matrix textures. We specifically discuss the stability of neutrino masses and lepton 
mixing angles with respect to the choice and variations of these small entries. 
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1 Introduction 



Recent data from the Super-Kamiokande experiment [|T], ||, |^ give a strong support to the 
interpretation of the atmospheric ||^ and solar ||^ neutrino problems in terms of neutrino 
oscillations. The results of the solar and atmospheric neutrino experiments along with the 
constraints from the reactor and accelerator experiments provided us with an important 
information on neutrino masses and mixing. This has resulted in an increased theoretical 
activity in the field with a large number of publications dedicated to the reconstruction of 
the neutrino masses and mixing and building theoretical models capable of producing the 
requisite neutrino mass matrices (see e.g. P, ||, ||, |lOi). 

In most of these studies some special structures or textures are suggested for the neutrino 
mass matrix. The special structures contain a number of large entries and also some small 
entries which are replaced by zeros in the textures. It is usually assumed that this form of 
the neutrino mass matrices emerge as a result of a flavour symmetry. In the exact symmetry 
limit some entries of the neutrino mass matrix vanish, but they are supposed to be filled 
in with small nonzero terms after a small breaking of the flavour symmetry is allowed. In 
many studies, however, these small entries of the neutrino mass matrices are not specified 
or are given only by an order of magnitude. 

While the neutrino mass textures reproduce correctly the gross features of the neutrino 
mass matrix, the small terms are of crucial importance for obtaining a realistic pattern of 
neutrino masses and mixing. This includes not only reproducing the correct values of the 
neutrino mass squared differences and mixing angles but also the correct correspondence 
between the ordering of the neutrino masses and the values of the mixing angles. Therefore 
any realistic neutrino mass model has to address the issue of the small entries of the neutrino 
mass matrices. 

In the present paper we consider the phenomenologically allowed forms of the neutrino 
mass matrices in the case of three light neutrino species, paying special attention to the small 
entries. The analysis is performed without resorting to any specific model of neutrino mass 
and therefore is quite general. Both hierarchical and quasi-degenerate neutrino mass cases 
are considered. Assuming maximal v^-Vr mixing, as suggested by the Super-Kamiokande 
atmospheric neutrino data, we derive simple approximate formulas which give the values 
of neutrino masses and remaining mixing angles in terms of the entries of the neutrino 
mass matrix. They allow one to find the neutrino masses and mixings off hand, without 
performing an explicit diagonalization of the neutrino mass matrix in each particular case. 
We discuss the stability of neutrino masses and mixing with respect to the choice and 
variations of the small entries of the neutrino mass matrices. The question of stability is an 
important issue in discussions of the radiative corrections in various neutrino mass models 
p. In the present paper it is studied purely phenomenologically, and the results can be 
used in the analyses of any models with three light neutrinos species. 
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2 Experimental input 



Currently there are three kinds of experimental indications of nonzero neutrino mass. These 
come from the solar neutrino experiments atmospheric neutrino experiments |I| 



and the accelerator LSND data [|TT|. The LSND result is the only one that has not yet been 
independently confirmed and we therefore choose not to consider it here. The results of the 
solar and atmospheric neutrino experiments can then be described in terms of only three 
light neutrino species, z/e, z/^ and z/^, without introducing sterile neutrinos. 

The solar neutrino deficit can be explained through the Ve ~^ oscillations where 
Vx is i^^i i'ti or their linear combination. There are three main types of solutions of the 
solar neutrino problem: small mixing angle MSW effect (SMA), large mixing angle MSW 
effect (LMA) and vacuum neutrino oscillations (VO). The atmospheric neutrino data can 
be explained through the oscillations with possible z/^ ^ Vf, oscillations as a 

subdominant channel. The corresponding allowed values of the neutrino parameters can be 
summarized as follows 1121 H: 
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(1) 

We shall use the parametrization of the unitary 3x3 matrix V describing the lepton mixing 
which coincides with the standard quark mixing matrix parametrization |1^. One can then 
identify the lepton mixing angles in ([|) as = 9i2, 6 atm = ^23- 

In addition to the results listed above, there is an important limit on the element t/es 



of the lepton mixing matrix coming from the CHOOZ reactor neutrino experiment 
|f/e3p(l — If^sP) < 0.045 — 0.02 for Am|^ = Am^j^ ranging in the Super-Kamiokande 
allowed region (2 — 6) ■ 10^'^ eV^, which together with the solar neutrino observations means 

sin^ ^13 = |f/e3p < (0.047 - 0.02) for Am^^ = (2 - 6) ■ lO^^ eV^ . (2) 

In the case of quasi-degenerate neutrinos we shall be also using the limit on the effective 
Majorana mass term mg// of z/g which comes from the Heidelberg- Moscow double beta decay 
experiment |]T5| : 

mejj<0.2eV. (3) 

This effective mass is essentially the irtee entry of the neutrino mass matrix. 

The results of the existing neutrino experiments can therefore be summarized as follows. 
There are two distinct mass squared difference scales, Amg and Am^^^, with the hierarchy 
AttIq -C Am^^^ ~ 10"'^ eV^. The mixing angle 6'23 responsible for the dominant channel of 
the atmospheric neutrino oscillations is close to the maximal one (45°), whereas the mixing 
angle 9i2 which governs the solar neutrino oscillations can be either small or large; for the 
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VO solution it can also be very close or equal to the maximal one. The mixing angle 
^13 responsible for the long-baseline Vf. oscillations and the subdominant v^^ ^ z/g 

oscillations of atmospheric neutrinos is either small or zero. 



3 Neutrino mass matrices 

The experimental information on the neutrino masses and lepton mixing angles allows one 
to reconstruct the phenomenologically allowed form of the neutrino mass matrix. We shall 
adopt the following convention. The neutrino mass eigenstate separated from the other two 
by the large Am^^^ is z/3 whereas those responsible for the solar neutrino problem which are 
separated by the small Amg are vx and v^. In other words, Arn^^^ ~ Am^^ ~ Amgg, 
Ami, = Am^i. Neutrino mass spectrum can be either hierarchical (mi,m2 <^ 1713 or 
mi, 7712 ^ rn^) or quasi-degenerate, mi ^ 1712 ~ m^, with only mass squared differences 
being hierarchical. The ordering of the z/i and 1^2 states is unimportant in the case of VO 
solution of the solar neutrino problem but it does make a difference in the case of the SMA 
and LMA solutions: the lower- mass state must have a larger z/g component. This condition 
will prove to be very important in deriving the constraints on the entries of the neutrino 
mass matrices. We shall disregard possible CP violation effects in the leptonic sector and 
assume the neutrino mass matrix to be real. Its eigenvalues, which we hereafter denote mj 
{i = 1,2,3), can be of either sign, depending of the relative CP parities of neutrinos. The 
physical neutrino masses are |mj|. 

We shall be assuming that the mixing angle responsible for the atmospheric neutrino 
oscillations 623 = 45° which is the best fit value of the Super-Kamiokande data f^. The 
case when 623 is close 45° but not exactly equal to this value can be treated similarly. In 
the first order in the small sin ^13 = e the lepton mixing matrix takes the form 
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i_ 

V 72( 



s + ce 
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6 \ 

1 

1 
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(4) 



where c = C12, s = S12, and the neutrino flavour basis is 



The neutrino mass 



matrix in the basis where the charged lepton mass matrix has been diagonalized can be 
written as tjil = t/diag(mi, m2, m-s) U'^ which gives 



rriL 



( 



73 [^("^3 



^[e(m3 - /i) + m_cs] ^[e(m3 - yu) - m_cs] \ 



/i) + m_cs] ^(m3 + fi' — 2m_cse) 
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(5) 



\{m3-ii') i(m3 + /i' + 2m_cse) j 



Here 



2 I 2 
mic -|- m2S , 



2 I 2 
niiS + m2C , 



m_ = m2 — mi 



(6) 
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We shall be using eq. (J^) to derive the phenomenologically allowed forms of the neutrino 
mass matrix for various neutrino mass hierarches. 

By making use of eqs. (||) and which are valid up to the small corrections of the 
order < 0.03, one can readily derive the neutrino mass eigenvalues and lepton mixing 
angles in terms of the entries rriij of the neutrino mass matrix rriL'- 



mi 2 



^ I ("^22 + "^33)/2 + mil - ± V [(^22 + m2,z) /2 - mil - "^23]^ + 2(mi2 - rnizY^ 

(7) 

mg = -(m22 + + 2m23) , (8) 

, \/2(mi2-mi3) A/2(mi2 + mi3) 

tan26'i2 = -tt ^ -, r , e = . (9) 

2(m22 + ^33) - (mil + "^23) m22 + ^33 + 2m23 - 2mii 

We shall now consider various cases of interest. 



3.1 Hierarchy |mi|,|m2| <C Ims] 

In this case the neutrino mass matrix (|^) takes the form 



rriL = mo 



e e 
e 1 + 5-5' 1-5 
\e' l~5 1 + 5 + 5' 



(10) 



with small k, e, e' , 5 and 5' . One can now use eqs. (0)-(0) to find neutrino masses and 
mixing angles up to the corrections of the order e^. It is not difficult, however, to find more 
accurate expressions, which include terms O(e^), directly from the matrix ( |T0|) 0: 
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13 



2^2(1 - k/2) 



(12) 
(13) 



^Up to possible trivial sign changes due to the rephasing of the neutrino fields. 

^ While (H) is valid only up to and including the terms 0{e), the form of the mass matrix (|l^) is generic 
for the case \mi\, \m2\ <^ \m3\, \e\ ^ 1. 
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where the eigenvalues rrii are given in units of tuq. The value of mo can be fixed through 
Am|^ ~ (2mo)^ — Am^^^. The mixing angle ^13 is small by construction and the limit 
d^) can be easily satisfied. The parameters Am2i = AmQ and du are fully determined 
by the small entries of the neutrino mass matrix, the values of which should be chosen in 
accordance with the assumed solution of the solar neutrino problem. 

For the MSW solutions (both SMA and LMA), the lower-mass state out of the two 
eigenstates, z/i and 1^2, must have the larger Ug component. The condition for this is 

\46-6'^\>\2K-{e^ + e'^)\. (14) 

Therefore the points in the parameter space at which the inequality in eq. ( [I^ ) is replaced 
by the equality are unstable: a small change of the values of the small parameters can affect 
drastically the pattern of neutrino masses and mixing. One of these points, 46 — 5'^ = 
2k — (e^ + e'^), is the point at which the denominator in the expression for tan 2^12 in 
eq. ( p!3D develops a pole, i.e. the mixing becomes maximal. Although at the other point, 
4(5 — 5'^ = — [2/t — (e^ + e'^)], the value of tan 2^12 remains finite, it is an instability point 
as well. This can be explained as follows. At this point one has m2 = —mi so that the 
eigenstates and V2 are degenerate. Crossing this point in the parameter space would 
interchange the relative position of the two eigenstates - the lower mass eigenstate would 
become the higher mass one and vice versa. At the same time, their flavour composition 
remains nearly constant since both 6'i3 and ^12 are regular at this instability point. This 
means that the condition that the lower-mass state out of the two eigenstates have the larger 
Vf. component is violated on one of the sides of this instability point. 

It should be noted that the expression 45 — 5'^ on left hand side of eq. (|14D is the 
determinant of the 2x2 matrix in the 23 subspace of the matrix (pUj). Therefore in models 
in which this subdeterminant exactly vanishes the condition ([T^ ) is not satisfied and the 
MSW solutions of the solar neutrino problem are not possible. 



3.2 Inverted hierarchy, |m3| <C |mi| 



(15) 



In the case jmsl -C |mi|, |m2| the neutrino mass matrix (|^) takes the form 

(KB e' \ 

e 1 + 5 — 5' —1 + 5 , 
e' -1 + 5 1 + 5 + 5' ] 

and the requirement |mi| ^ |m2| leads to some additional constraints. There are essentially 
three possibilities: 

(1) e, e' 5 and 5' small, k ~ 2. This gives mi ~ m2. The neutrino mass eigenvalues and 
mixing angles are given by 



mi 2 ^ mo 



1+-±W 1-- + 



[e — e 



/^2 



ms 



25 - — \mQ 



(16) 



5 



(e -e') e + e' 

tan2^i2~^:^ — , sin^i3 = e - ^ . (17) 

v^(l - k/2) 2y/2{6 - Kjl) ^ ' 

(2) £, e' 8 and 8' small, k ^ —2. This gives m2 — —mi. Formulas (p!6[) and (p!7[) apply 
in this case, too. 

The case 1712 — —mi obtains also when e and e' are not small provided they are nearly 
negative of each other. In this case {e' ~ —e, \e\ ^ 1, 6 and 6' small, k, ~ —2) the values 
of mi^2 and tan26'i2 are again given by the first equations in (|T6|) and (p!7|), whereas for 
and sin 613 one obtains 

{e + e'y + {e- e'^S + {e^ - e'^)S' - AnS 

(.^ + .-^)- 2.(1 + ^) -4^ ^''^ 

sm^i3 . is + e')-6ie-e')-e'i6' + m,/m,) ^^^^ 

Je'^{e'^ + e'2) — 2£:(£ — e'){K — ms/rrio) + 2(/t — m^/mQY 



atm ' 



In both cases (1) and (2) the value of tuq can be fixed through m^ ~ mg — (2mo)^ — Am; 
whereas Am2i = Amg ~ 8mgr where r is the square root that appears in the first equation 
in (plGl). The mixing angle 613 is small by construction; in case (1) the value of tan26'i2 
depends on how close k is to 2 as compared to e — e'; in case (2) tan 26^12 is small when e 
and e' are small but can be large in the case e' ~ —e, \e\ > 1. 

The condition that the lower-mass state out of the eigenstates ui and 1^2 have the larger 
z/g component, which is important for the SMA and LMA solutions of the solar neutrino 
problem, is 

\k\<2. (20) 

It constrains only a large entry of the mass matrix and not the small entries. This is related 
to the particular choice of the parametrization of the mass matrix in eq. ([T5| ) which is 



convenient because it leads to the very simple expressions (T^) and (P!7|). In general, this 
condition reads 2|mii| < |m22 + mss — 2m23|, which includes both small and large entries. 
The parameter k is required to be close to ±2 in the cases under discussion. A small change 
of its value can drastically alter the pattern of the neutrino masses and mixing. As an 
example, consider the case k = 1.998, e = 5 - 10"^, e' = 1- 10"^, 6 = 1- 10~^ and 5 = 3- 10~^. 
This yields mi ~ 1.996, m2 ~ 2.002, rris ~ 0.02 (in units of mo), sin 6*13 ~ —2.14 ■ 10~^, 
cos6'i2 — 0.816 (sin^ 2^12 — 0.889). This choice of parameters is suitable for the LMA 
solution of the solar neutrino problem. The fact that cos^ 6'i2 > 1/2 means that the lower- 
mass state z/i has the larger z/g component than the higher mass state z/2, as it should. If, 
however, the value of k is increased by just 0.2% while the other parameters are kept intact, 
the eigenstates z/i and 1^2 interchange their position so that now mi > m2 (alternatively, 
one can rename the states in such a way that mi is always smaller than m2 but then sin ^12 
and cos ^12 interchange). This means that the lower-lying of the two mass eigenstates no 
longer has the larger Uf. component and therefore there is no solution to the solar neutrino 
problem at all. Similar situation takes place in case (2) when n ~ —2. 
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We have seen that, hke in the case of the normal hierarchy |mi^2| ^ I'^al, in cases (1) 
and (2) of the inverted hierarchy there are instability regions which are close to the certain 
points in the parameter space. In the case of the normal hierarchy these instability regions 
are just small domains of the full allowed parameter space, and the parameters can in general 
be quite far from these domains. On the contrary, in the cases of the inverted hierarchy 
under discussion, the parameter n must be close to ±2, so that the situation is generically 
unstable. This is related to the fact that the states z/i and 1^2 are quasi-degenerate in this 
case while in the case of the normal hierarchy their masses may be hierarchical. 

(3) e' ~ ±e, \e\ ^ 1, \6\, \6'\; \k,\ 1. 5 and 6' need not be small. This case also leads to 
m2 — —nil. By a rescaling of the parameters it can also be formulated as e' ~ ±e = 0{1), 
the rest of the entries of the matrix being small. This case was suggested in on the 
basis of the approximate — L^j_ — symmetry. It is more convenient to use a different 
parametrization of the mass matrix in this case: 

j e 1 ±1 \ 
mL = mo 1 5 e' , (21) 
V ±1 e' 6' J 

where e, e', S and S' are small. The neutrino masses and mixing angles are now 

l^^[2{e ± e') + 6 + 6'] - ~ mo |i[2(£ ± e') + ^ + 5'] + v^j , (22) 

m3^^(5 + 5'T2e'), sin^ia^^T^- (23) 

where the upper and lower signs refer to the corresponding signs of the matrix element 
(m-L)i3 in (|2T|). In this case, too, the mixing angle 613 is automatically small. The mixing 
angle 612 = 45° up to corrections of the second order in the small entries of the mass matrix, 
i.e. only the VO solution of the solar neutrino problem is possible 0. This is an example 
of the so-called bimaximal mixing p|. Unlike cases (1) and (2) discussed above, this case 
is very stable with respect to variations of the entries of the neutrino mass matrix. This is 
because the maximal mixing 612 = 45° makes the ordering of ui and z/2 unimportant. 



mi ~ mo 



3.3 Quasi-degenerate neutrinos, |mi| ~ \m2\ — \ms\ 

The quasi-degenerate cases can be classified according to the relative signs of the neutrino 
mass eigenvalues. 

(1) mi ~ 7712 — m3. We shall define the mass splittings 6 and A through m2 = mi -f- 5, 
m3 = mi + A, \6\ ^ |A| ^ |mi|. Then the neutrino mass matrix (^ can be written as 

niL = midiag{l, 1, I) + , (24) 

^This remains true also in a more general case in which (771^)13 — A, (mi)i3 — B and for A ^ B the 
mixing angle ^23 — ~ arctan(i3 /A) is large but not necessarily maximal. 
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where the matrix M has the same form as the one on the right hand side of eq. ([ToD . 
Therefore the formulas of sec. 3.1 apply to this case up to the obvious modification of the 
mass eigenvalues. Zeroth order texture in this case is proportional to the unit matrix. 

(2) mi ~ 1712 — ""^3- We now denote m2 = mi + 5, = —mi + A. The mass matrix 
(I5|) takes the form 



rriL = mi 



/ 



V 



1 -V2e 
-V2e 
-V2e -1 




+ 



5 
2 



2cse 

2 



(25) 



2s^ 

V2{cs - es2) 

-V2{cs + es^) -c^ + 2cse 

where the notation is the same as in eq. (^). The zeroth order texture is given by the first 
term in this equation in which e is set equal to zero. Neutrino masses and mixing angles 
can be expressed through the entries of the mass matrix (|25|) with the use of eqs. (0)"(il)- 

(3) mi ~ —1712 — —TTT'3 0- We now denote m2 = —mi + 6, = —mi + A. The mass 
matrix (^) takes the form 

mi = Ml + Ma + Ms , (26) 

where M2 and M3 coincide with the second and the third terms on the right hand side of 
eq. (P3|) respectively, and 



/ 



Ml 



-mi 



V2{cs + ec2) 



V2{cs + ec^) -V2(cs-ec2 



2cse 



(27) 



V -V2{ 



cs — ec 



& + 2cse 



Notice that, unlike in the two previous cases, the zeroth order texture (which obtains from 
([27| ) in the limit e = 0) depends on the mixing angle 6*12. This is because the exact degeneracy 
limit now corresponds to m2 = —mi and not to m2 = mi, so that the mixing between Vx 
and V2 remains meaningful. Therefore the zeroth order texture in this case depends on the 
choice of the solution of the solar neutrino problem. For the SMA solution one has to choose 
c = 1, s = or c = 0, s = 1; the corresponding zeroth order textures are diag{—l, 1, 1) and 
the matrix with mn = m23 = m32, the rest of the elements being zero. Small nonvanishing 
^12 results when the zeroth order textures are perturbed. Two further examples are 



/ 
V2 
V V2 



V2 
1 
1 




1 
3 



-1 
2 
-2 



:28i 



the first of which corresponds to the bimaximal mixing (^12 
an example of the texture with a large but not maximal 6*12 

^The case mi ~ — m2 



= 6*23 = 45°) and the second is 
s = 1/^3 (sin^ 2^12 = 8/9). 

TO3 can be obtained from the present one by the change of notation vi ^ 1/2- 
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We shall now consider one of the examples, namely, c ~ 0, s ~ 1, in more detail. The 
neutrino mass matrix for this case can be written as 



mL = mA e 6 1 + 6' \ (29) 




with small e, e', S, S' and p. From eqs. (0)-(|^) one finds 



(30) 

e — e' e + e' 

tan 2^12 = 1= , sin0i3 = —pz . (31) 

V2[l + 5V2-(p + 5)/4] ' y2[(p + 5)/2 + 5'] ^ ' 

The mixing angle ^13 depends on a ratio of the small parameters. It is not automatically 
small and care should be taken in order to satisfy the constraint (0). The mixing angle 612 
is generically small because the zeroth order texture corresponds to zero mixing. Thus, in 
this case only the SMA solution of the solar neutrino problem is possible. The neutrino 
mass splittings are 

5 ~ mo [(p + 6)/2 -6'], A ~ mo [p + 5] . (32) 

The condition that the lower-mass state out of the eigenstates z/i and 1^2 have the larger 
component is 

5'>ip + 6)/2. (33) 

The point at which the inequality is replaced by the equality is unstable with respect to 
small variations of the parameters. As follows from (^2]), the requirement of the correct 
hierarchy of the mass differences, |A| ^ \6\, is equivalent to p + 5 ~ 26'. This means that 
in the present case the parameters are necessarily close to the instability point. 

The value of the mixing angle Oi^ as well as those of the neutrino mass splittings depend 
sensitively on the value of S'. In order for 6/mo and A/mo to be small so must be 5', which 
we assume. Notice, however, that 6' is the difference of the two large matrix elements m23 
and mil. This means that a very small relative change of these elements can result in a 
drastic change of the whole pattern of neutrino masses and mixing. This situation is typical 
for all cases of quasi-degenerate neutrinos. 

Finally, it should be noticed that if neutrinos are quasi-degenerate, care should be taken 
in order not to violate the limit (|^) coming from the experiment on the neutrinoless double 
beta decay. It essentially constrains the mn entry of the mass matrices. From the discussion 
in this subsection it follows that the common mass scale mo must satisfy this constraint in 
all quasi-degenerate cases except in the one of bimaximal mixing, for which the zeroth order 



texture is given by the first matrix in eq. (28). This case leads to the VO solution of the 
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solar neutrino problem. As follows from the second example in (^), the bound on mo can 
also be somewhat relaxed in the case of the LMA solution of the solar neutrino problem. 
In all quasi-degenerate cases in which the common mass scale has to satisfy the limit (|^) 
neutrinos cannot constitute any appreciable part of the dark matter of the universe. 



4 Conclusion 

We have studied systematically the phenomenologically allowed structures of the neutrino 
mass matrix in the case of three light neutrino species. For both hierarchical and quasi- 
degenerate cases we considered neutrino mass matrices which can be obtained by perturbing 
the well known zeroth order textures. 

Assuming maximal i'^-i't mixing we derived simple analytic formulas which give neutrino 
masses and the remaining mixing angles in terms of the entries of the mass matrices. We have 
checked our approximate analytic expressions by performing the numerical diagonalization 
of neutrino mass matrix and found a very good agreement in each case. 

We analyzed the stability of the neutrino masses and mixing angles with respect to 
small variations of the entries of the mass matrices, paying special attention to the small 
entries. In particular, we have studied the stability of the condition of having a larger 
component in the lower-mass state out of the eigenstates ui and 1/2 responsible for the solar 
neutrino oscillations. This condition is important for the SMA and LMA solutions of the 
solar neutrino problem. 

We have found that while the normal hierarchy case \mi\, |m2| <C jmsl is in general 
very stable, the inverted hierarchy case Im^l <^ |mi| ~ \m2\ is much less stable (except in 
the case of bimaximal mixing discussed in case (3) of sec. 3.2). The quasi-degenerate case 
suffers from a number of instabilities and therefore is the least natural one. 

The author is grateful to G. C. Branco and M. N. Rebelo for useful discussions. This 
work was supported by Fundagao para a Ciencia e a Tecnologia through the grant PRAXIS 
XXI/BCC/16414/98 and also in part by the TMR network grant ERBFMRX-CT960090 of 
the European Union. 
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